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An alternative approach to the derivation of the force on a vortex based in an adiabatic approx- 
imation in the action of the superfluid system is developed. Assuming that the vortex motion is 
relatively slow compared with the characteristic times involved in the microscopic degrees of freedom, 
the effect of the superfluid and its excitations is reduced to a gauge potential, and the associated 
transverse force. The vortex velocity part of the transverse force is found in terms of the thermal 
expectation of the angular momentum of the fluid around the vortex. The excitations countercir- 
culate the vortex reducing the effective density to that of the superfluid part only. Non-adiabatic 
\^ , contributions appear in this formalism as a non-Abelian gauge potential connecting different mi- 

^\ • croscopic states, in particular the renormalization of the vortex mass is given by the lowest order 

0^ I diagonal correction to the adiabatic theory, and the long wavelength contributions found to depend 

only on the static structure factor for the fluid. 
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INTRODUCTION 



> . . 

■ Quantized vortices have been an essential part in the theory of superfluids since Onsager first stated the idea of 
' quantized circulation in the late 40'a2l'a. In fact, vortices in general have been an important part of classical fluid 

mechanics for decadesd, in particular in the study of turbulent regimes. The Magnus force , or Kutta-Joukowski 
hydrodynamic lift, Fm = p k x {vy — I'fluid) is well knnwn from classical hydrodynamics and occurs whenever an 
object with circulation k around it moves through a fluidQ, an important application being the lift force on an airplane 
wing. The coupling to the fluid is also known to produce a hydrodynamic mass. 

The problem of obtaining effective dynamics for adiabatic and non-adiabatic motion of vortices is highly controver- 
"j^ • sial. In particular, the expression for the Magnus force at|-finite temperatures is still far from clear, and today (more 

■ than three decades-after the first measurements by Vinerfl), different expressions for the Magnus force can be found 
I 1 in the literatureUllj. The confusion arises, in part, from different interpretations on the role played by excitations 

''O being scattered asymmetrically at the vortex, leading to a transverse force proporticpial to the normal fluid density 
Ch , Pn and either the relative velocity (t>n — vy) or {v,^ — Vs), namely the lordanskii forcalj. The magnitude of this term 
O ■ must be calculated and, moreover, it is not clear whether it should be added to the Magnus force written above, with 
' the coefficient p, or to a similar expression involving the superfluid part ps only. This paper analyzes the situation for 
^ ' neutral, homogeneous superfluids, while the case of charged superfluids and non translation invariant systems will be 
'k>( [ published elsewherell3. ■— ■ P-, 

For fermion superfluids, VolovikEJ and StoiieEJ have identified additional contributions from localized quasiparticles 
inside the vortex core. Thouless, An and NiullJ(TAN) show that no such a contribution exists for translation invariant 
systems, as seems to be observedll3 for the B phase of ^He. 

Another point of controversy is the effective mass. Coupling of the vortex motion and the superfluid circulating 
around it renormalize the masS|-and it has been argued by some authorSpUaat this renormalized mass is logarithmically 
divergent with the system sizalJ, while others estimate it to be finitdlS or even zeroll3. Recent studies by us and 
Demircan et aLt3 show that the vortex mass is logarithmically dependent on the frequency, cutting off the size 
dependence for large systems. 

In section ^ we develop a formalism for the effective action of a vortex moving relative to the fluid. Following TAN, 
this is accomplished by the introduction of a short range repulsive potential for the particles, strong enough to pin 
the vortex (a simple example would be a macroscopic solid wire as doiie in Vinen's experimentO) . We follow closely 
the adiabatic approximation as done by Moody, Shapere and WilczekEll, and arrive to a simple exact expression for 
the effective action involving the vortex coordinates. Expansion of this action yields the vortex velocity part of the 
Magnus force (VVPMF) to first order in the vortex velocity (sec. ||). This force depends on the total circulation 



around the vortex, which is calculated in section [IT The excitations countercirculate the vortex in a way that is 
equivalent to stating that only the superfluid fraction has a non- vanishing circulation as stated in TAN. The VVPMF 
is given by psK^ x vy- 

The normal and superfluid velocity parts of the Magnus force (NVPMF, SVPMF) require extra work, and will be 
dealt with in a future publication. 
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Section IV deals with the second order terms, which renormahze the vortex mass. This mass renormahzation 
depends only on the static structure factors for the fluid. We recover a logarithmically divergent mass in the zero- 
frequency limit. 



I. ADIABATIC EFFECTIVE ACTION FOR A VORTEX 



Following reference TANt3, we control the motion of a vortex by means of a pinning potential, consisting in a 
short-range repulsive potential for the particles that form the fluid. At this stage we only require this interaction to 
be a pure potential V{R,r]) without any velocity dependent part, where R iSj-the vortex coordinate and rj = {r,} 
is the set of all particle coordinates. We follow Moody, Shapere and WilczckO rather closely, avoiding some sign 
inconsistencies along the way. We can immediately write the Lagrangian of the system as 

L = ibarc + 'fluid(-R, 77), (1) 

where Lbarc is the Lagrangian of the bare pinning center and lnuid incorporates the full Lagrangian of the fluid and 
the repulsive potential V{R,ri). Notice that we have imposed no condition on the form of the Lagrangians, except 
for the type of interaction between the pinning center and the fluid particles. 

The full tirneievolution kernel of the system can now be written in terms of Fcynman path integrals over all possible 



configuration; 



J22 



U{Rf,7]f,tf;R,,r],,U) = {Rf,Tjf,tf\R^,T],,U)=J^J P[7?] exp |- j [U,,, + huidiR, v)]dt j ■ (2) 

For a stationary vortex at position R we can also, in principle, obtain the exact eigenstates of the fluid subsystem 

hiR)MR) = ^n{R)MR), (3) 

and for any particular value of R, any wavefunction for the fast system can be expanded in this complete basis 

^ = ^(/.„(H)^^„(H). (4) 

n 

We will be interested in finding a kernel Umn{R}, tf', Ri, ti) that links the final weights Fm{Rf, tf) with the initial 
state's Fn{Ri, ti): 

Fm{Rf,tf) =^t/™„(i2/,t/;i2^,i^) Fn{R^,ti). (5) 

n 

This kernel has all the information that we need, and is diagonal in the stationary case R = 0. Most important is 
the fact that, for non-degenerate cases, off-diagonal contributions are small. This new kernel is given by 

Umn{Rf,tf;R^,ti)= jj drjf dr]^ {4>,ntf\-nf) U{Rf,rif,tf;Ri,ri^,ti) {rii\4>nti) 

^'v[R] K'^-'-U^^", (6) 



(7) 



iR^ 



Vi 



We can calculate this propagator in the usual way by dividing the time interval [ti, tf] into N segments and inserting 
the identity operator at each time. Taking N ^ 00 yields 

I, I., , ,, — 1 v_i ■' 



Iu...,In-i3 = 1,---,N-1 

Texp I J f ' dt [~e,n Smn+ifl R{t) ■ {(f>,n\'V B.<f>n)] } , (8) 
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where the last expression is formal way of denoting an infinite number of time ordered matrix products. 
The time evolution of the pinning center can now be written in terms of an effective action 



rRf 

Umn{Rf,tf;R,,U)= ' V[R]T 6^^°"-, (9) 

J Ri 



rtcff _ 



dt [Lb 



+ in R- 



(10) 



The results insofar are exact, but not very helpful. In general, the gauge-like term ihR - {4>m\^ Rf^n) is non-Abelian, 
and we have to deal with complicated time ordered products, etc. For non-degenerate states, however, it is possible 
to show that non-diagonal terms vanish faster than any power of {TAe/h)~^ , with t being a characteristic time of 
motion of the pinning center, and Ae the smallest energy difference of the fluid's dynamics (see ref. 21). In most cases 
it will be possible to greatly simplify the problem by only dealing with a few degenerate or closely degenerate states. 

For non-degenerate cases, transitions to other states may be neglected and we can easily calculate the effective 
action for a given state of the fluid to any order in the vortex velocity. In particular, the Magnus force is linear in the 
velocities, and effective masses correspond to quadratic terms in the action. Expanding to these orders, the effective 
action for the vortex motion is simply given by 



dt 



(11) 



{di(j)n\(f>l){(l)l\dj(j)„ 



(12) 



II. VORTEX VELOCITY PART OF THE MAGNUS FORCE 



From the effective action (|11| ) , we can immediately find the transverse VVPMF for any given quantum state 0„ of 
the fluid: 



ihkx {Vr X ((^„|Vj^0„)) . 



(13) 



Identifying vy — R, we now perform a statistical averaging over initial states to obtain the thermally averaged 
VVPMF 



F X z= -i h Vv fn 

n 

= -ihvv ^ /„ 



dX' dY 



dY ' dX 



(14) 



where /„ is the probability of occupation of the state This is the familiar form of the Berry curvaturaHj. For a 
homogeneous system it is possible to substitute V/? = — Vi, and therefore 



■ '-'Vn I '->Vn ^ 

dxi duj ' 



dm dxj 



ihz-J d\[V xV'p{r\r)]r^r' - iTij d^r j rfV[2Vi x V2'r(ri', ra'; n, ra)]. 



(15) 



The last term in the equation above vanishes, since the fitst line above is the commutator of the x and y components 
of the total momentum which is a one particle operatoi-Eij. The integrand of the first term is equal to V x (V — 
'V')p{r,r')/2. Application of Stokes' theorem can be used to write the force per unit length 
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F/vv = - y j( di-[(V - V')p{ry)]r=r' 



(16) 



This result is exact, and the contour of integration T can be taken as far from the vortex as one wishes, and there 
are no contributions from the vicinity of the core. 

For a neutral superfluid the integrand is just the momentum density j. At this stage one can decompose the 
momentum density in a two fluid picture as j = psVg + Pn I'n, and the vortex velocity part of the Magnus force will 
be given by the superfluid and normal mass densities times the superfluid and normal circulations 

F/vv = <j> {PsVs + PnVn)-dl = Ps Kg + Pn Kn- (17) 

The circulation of the superfluid Kg is quantized to multiples of h/m, and both ps and p„ are well defined quantities. 
TAN have argued for the normal fluid to have no circulation. However, the circulation of a normal fluid is something 
that has to be explicitly calculatedEa, which is done in the following section. 



III. CIRCULATION OF THE NORMAL FLUID 



As stated before, one cannot in general assume the value of the circulation of a regular fluid. Aeronautical engineers 
put up a fair amount of effort to calculate the circulation around a particular wing section.|-JFortunately our case is 
considerably simpler. We follow the method used by Landau to calculate the normal densitynS in order to obtain the 
circulation of the normal fluid. The VVPMF is then completely defined in terms of macroscopic transport parameters. 

Note that the expression for the Magnus force is given in terms of the circulation, and that it can be defined by any 
path circling the vortex. Consider a vortex centered in a cylindrical container and "average" over all possible paths 
that go once around the vortex 

F/vv [ dz^ [ rdr (f j-dl = 27r — ^ d^'x z-{r x j) = 2n^. (18) 

Lz Jo R Jo Jr(z,r) ttR^Lz V 

The Magnus force per unit length is just 27r times the average angular momentum Cz per unit volume V. At zero 
temperature, all the fluid is superfluid, each particle averages and angular momentum fi. At zero temperature (or in 
the ground-state), the Magnus force is simply given by 

UN h 

FoW = 2tt — ^ p -. 19 
V m 

Finite temperatures can be analyzed by adding the effect of excitations. At low temperatures excitations can be 
dealt in a dilute gas approximation. Moreover, phonons become predominant and we may neglect rotons. At this stage 
one wishes to calculate the expectation value of the angular momentum of the excitations. The presence of the vortex 
couples the excitations to the superfluid velocity field, Doppler shifting the energy of excitations, thus increasing the 
occupation of countercirculating modes relative to those that circulate in the direction of the superfluid. 

The expectation of the angular momentum of the excitations is simply given by thermally populating the different 
phonon modes (labeled in cylindrical coordinates by the radial and axial wave- numbers, and by the azimuthal angular 
momentum): 

The mode frequencies can be obtained from wave equations derived from either a Feynman many-body stateEl or a 
non-linear Gross-Pitaevskii theorycH. From this latter point of view, we start with a non-linear Schrodinger equation 
to describe a weakly interacting Bose gas, and after some algebra (Appendix ^) get the following wave equation for 
phonons in presence of a static vortex: 

c2 - V' - — <^ • V^A = 0. (21) 
Trr 

This last term decreases the frequency of phonons countercirculating the vortex. In a cylindrical geometry, the 
eigenstates can be labeled by the angular momentum, radial and azimuthal wave- vectors. To a good approximation 
the normal mode frequencies are given by: 
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C 



.v/fefTIf + m-^ (-). (22) 



27r V2' 

Using this information in equation (^0|) we get to first order in the superfluid circulation 

where now the sum is performed on the normal modes in absence of the vortex. We now assume that the typical 
inter-phonon equilibratiotuidistance is much smaller than the size of our container, which might be rather unrealistic 
at very low temperaturesEa, but makes sense for a theoretical infinite volume, infinite time-scale limit. In this limit, 
equilibration occurs locally, rather than at the boundaries of the system, and we may simplify the summation over 
modes by a sum over uniformly distributed states in phase space {J2k k m ^ I ^'^^ '^^^ /(^ti")"^, m ^ (r x fe)^ and 
< l/r^ >^ l/r2): 

ft2 n27T poo j-oc hck/ksT o 2 (U^rp^^i 



" ('^^^^^ i I ier^^^/^sT _ 1)2 = -(-^^^^) (24) 

which shows the excitations are indeed giving some countercirculation around the vortex. The total angular momentum 
is just given by the sum of the ground state (or zero-temperature) angular momentum and that of the excitations, 
which for Kg = h/m yields 



_ ttR^L^ h 

•C — £0 + -Ccxcit — n — 

Ztt m 



27r2 (kBTy 



45 n'^c^ 



nR^L, h nR^L, h 

— [P-Pn\ = ^ Ps, 25) 

ZTT m ZTT m 



where the last two identities follow from identifying the second term in the brackets as the normal fluid density due 
to phonona2j. This magnitude of the excitation countercirculation is equivalent to the non-circulation of the normal 
fluid in the two fluid model. Application of equation (|l^) gives the correct vortex velocity part of the Magnus force 

h , , 

-Fvortox ^ Ps — Z-KVy. (26) 

m 

Where the result above was calculated with — — 0. Next section deals with the coefficient of the force 
proportional to that is the NVPMF. 



IV. VORTEX MASS 



In this section we calculate the renormalization of the vortex, mass in the static limit. Vortex mass calculations 
have been performed starting from the hydrodynamic equationstj, and here we approach the problem from a different 
(albeit equivalent) point of view: from the second order correction to the adiabatic approximation in the effective 
action (eq. |ll| and ^2|). This second order adiabatic approximation yields correct results for the long wavelength 
contribution to the vortex mass. |-||— ■ 

For a superfluid, one can write Feynman's many-body wavefunctionacj for a static vortex at R as 

i 

where |(/)o) is the many-body ground state. In this simplified case 

|5^^) = E^^ !'/'>= \dv^^Y.^{T-T^)\i,)= [dr^^pirm, (28) 
xf +yf J x'^ ^ J x'^ + 

where p[r) is the density matrix. This can be Fourier transformed to read 

/dk' k' 
J^A^'^)^2Pk' 1^)- (29) 

It is clear that the only significant couplings in eq. (|l^) will be to states that differ from the "one vortex ground 
state" by p, that is the one vortex plus one phonon states: 
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|fc)«Pfe|V). (30) 

We can readily calculate the overlap 

/dk' k' 
^ 5{k'^)^{^\plpj^,\i^) I (^l^pfcl^)^/^ (31) 

Under the approximations used so far the correlation in the last term can be substituted by the ground state 
correlation 

{^\piPk'\^) « {MpiPk'\M=N S{k)5^j^,, (32) 

Where S{k) is the static structure factor!! and, in neutral superfiuids, is equal to fc/(2mc) in the long wavelength 
limit (c is the speed of sound). We are left with 



1/2 

(33) 



Given that e(fc) « /c^/(2m5(fc)), we can write 

= = 0, (34) 



Mr = Mr = ,n^pL. [ d^k 'Isikf = ^ /'^"'"^ f = ^ log (^) , (35) 

Where R and a are large and small distances cutoffs. In the static limit one obtains the size dependent logarithmic 
divergencejif the mass as predicted by DuarO from the hydrodynamic equations of the supcrfluid. A finite frequency 
calculations will cut R, by the phonon wavelength R^ \ — 27rc/a;. 

A very illustrative fact is that the additional vortex mass is related to the energy of a static vortex by the Einstein 
relation i^static = Mq . Where 

i^static ^ fdrlpv^^^ log (^) . (36) 
J 2 m'' V a / 

This is not surprising, and should be expected for any wave-equation like dynamics like the ones involved for 
phonons. 



V. CONCLUSIONS 



Starting from an adiabatic expansion for the effective action of a vortex we have obtained, in a closed and coher- 
ent way, the vortex velocity dependent part of the transverse force acting on a rectilinear vortex in a neutral and 
homogeneous superfluid. 

In this circumstances, the force that depends on the vortex velocity can be completely determined by conditions far 
away from the vortex core and is proportional to the superfluid density, which is correct in light of Vinen's experimentQ. 

Within the formalism used, the renormalization of the vortex mass depends only on the form of the static structure 
factor for the fluid. The obtained result coincides with similar calculations derived from hydrodynamic equations. 

More research is currently under way to completely determine normal and superfluid velocities dependent part of 
the transverse force, and to elucidate the effects caused by charge and the loss of translation invariancal3, as well as 
the roles of phonon radiative damping and the longitudinal part of the lordanskii force. 
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APPENDIX A: WAVE EQUATION FOR PHONONS IN PRESENCE OF A VORTEX 



E 



Starting from the Gross-Pitaevskii nonlinear Scrodinger equatior 

*^7?^-T-^'* + ^(l*l'-^o)*, (Al) 
ot 2m 

a Madelung transformation0 4* — ^Jp e^^ leaves us with the continuity equation and the Euler-BernouUi equation^: 

|f + V-(pv)=0, (A2) 



Is+^^ + c^P^^Q. (A3) 

Where v = fi/m and we have left out the "quantum pressure" term, which is irrelevant for our purposes. After 
some algebra, separating the static phase due to the vortex and linearizing, one obtains the following equation of 
motion for the phonon part of the phase: 

V^V - V' - —0 -V^/' = 0. (A4) 

This expression is not valid in the immediate vicinity of the vortex, where the assumptions made are not valid. In 
general this will only be important when either dealing with short-wavelength excitations or the s-wave modes, that 
can reach these regions. 
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